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Abstract. The ATSP polytope can be expressed by asymmetric polynomial 
size linear program. 



Introduction 

In his seminal work [1], M. Yannakakis proved that the Traveling Salesman 
problem (TSP) polytope cannot be expressed by symmetric polynomial size linear 
program, where symmetry means that the polytope is an invariant under vertex 
relabeling. The question about the size of asymmetric linear programs was left 
open in [1], and it remained open since. 

This article answers the question. We reduce the Asymmetric Traveling Salesman 
problem (ATSP) to an asymmetric polynomial size linear program whose feasible 
set is asymmetric in the sense of pQ. 

1. Directed Hamiltonian Cycle Problem 

Directed Hamiltonian Cycle Problem (DHC) is a problem to find in any given 
digraph a cycle which visits all vertices (a circuit which visits all vertices and visits 
each of them just once). Any such cycle is called Hamiltonian. 

Ultimately, the problem can be solved by testing all vertex permutations on 
being a Hamiltonian cycle, or by "growing" of paths. The only drawback here is 
the number of "elementary operations" required in worst casqj- Computational 
complexity of the methods will be factorial and exponential appropriately. So, they 
are infeasible for modern computers even in the cases of modest digraphs. 

The real problem is to detect Hamiltonian cycles in feasible time. Sometimes, 
that can be done using the appropriate criterion [H [31 HI [SJ [51 [71 [HI and many others]. 
But, such particular instances of DHC are rare. 

The importance of DHC follows from the fact that its decision version^ is a NP- 
complete problem^ [9]. There is a wide spread belief in P ^ NP. This hypothesis 
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x That number as a function of digraph's size is called a time or computational complexity. 
2 The problem of the existence of Hamiltonian cycles. 

3 NP-complete problem is such NP-problem to which any other NP-problem can be reduced in 
polynomial time by deterministic Turing machine (DTM). DTM is a rigorously defined algorithm. 
NP-problem is a problem which can be solved in polynomial time by non-deterministic Turing 
machine (NDTM). NDTM is an aggregate of a non-deterministic polynomial time generator of 
guesses for solution and a polynomial time DTM which checks those guesses on being a solution. 
P-problem is a problem which can be solved by DTM in polynomial time. For more information, 
see, for example, |10| . 
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means that no polynomial time algorithm exists for DHC [9j [10] . Yet, below we 
express DHC by a polynomial size linear system. It is well known that solution of 
such systems is a P-problem 

1.1. DHC as a relabeling problem. Let G be a given digraph with n vertices, 
ri > 1. Let's arbitrarily enumerate/label vertices of G. Let Aq be the appropriate 
adjacency matrix of G: 

Ag = \9ij)nXn 

- where gty is or 1 depending on the existence/absence of arc from vertex i into 
vertex j. 

Obviously, there is a Hamiltonian cycle in digraph G iff there is a circular per- 
mutation submatrix in matrix Aq. That submatrix is an adjacency matrix of that 
cycle. This trivial fact may be expressed by the following quadratic systerrQ: 



(1.1) 



Aq > XSX 1 



- where unknown X is a permutation matrix, and constant S is any circular per- 
mutation matrix - an adjacency matrix of a Hamiltonian cycle. Let's fix to S the 
following value: 
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Any permutation-matrix solution X of system 11.11 means a Hamiltonian cycle 
in G. Matrix XSX T is an adjacency matrix of that cycle, and matrix X T defines 
relabeling of G along it. After that relabeling, this cycle will visit vertices of G in 
the following order: 

1^2^3^...^n^l 

System [TTT1 presents DHC in all senses: DHC as a decision problem is the con- 
sistence problem for the system; DHC as a search problerrQ is the problem to solve 
the system; DHC as a counting problerr@ is the problem to count the system's so- 
lutions, etc. 

In terms of system 11.11 guesses are vertex permutations. They are presented 
in the system with unknown permutation matrix X. Each solution of the system 
delivers a Hamiltonian cycle. 



1.2. Compatibility matrix. In system ITTTl digraph G is presented with its ad- 
jacency matrix. Adjacency matrix is a universal encoding for digraphs. To solve 
system 11.11 we will use special encoding which we call a compatibility matrix. 



For two matrices B = (bij) and C = (cij) of the same size, relation B > C means that 
Vi,i (bij > Cij). 

^The problem of the actual finding of Hamiltonian cycles. 

^The problem to count different Hamiltonian cycles in G - n times of that number is the 
number of solutions of system 11.11 
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For each two vertices i and j, let's build a compatibility box. The box is the 



following matrbfl dj = (c yA1I/ )„ X n: 
(1.2) 



The (p., ^)-th element of Cy indicates whether vertex couple ([A, v) can be relabeled 
into couple (i,j) regardless of anything but adjacency. 

There are n 2 compatibility boxes, and it takes time 0(n 2 ) to compute any of 
them with brute force. The boxes have the following major properties: 



(1.3) 



^ j => Ci jflfl = 



- where U n is n x n identity matrix, and all indexes are in their ranges. 

Compatibility matrix is a box matrix which aggregates all compatibility boxes 
in accordance with their indexes: 

C = (Cjj )nxn 

In elements, the matrix has size n 2 x n 2 . 

Compatibility matrix aggregates all relabeling options for vertex couples of G. 
For our particular value of S, the matrix looks: 
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\ A G K n K n 
where the size of C is shown in boxes and 

Kn = ( 1 ) n. X 



A T 

G 



U n 



U n J 



Obviously, the pattern holds for any circular permutation matrix S. 

Due to the previous subsection, G has a Hamiltonian cycle iff its vertex couples 
can be relabeled without any contradiction, i.e. iff C has a grid of non-zero elements, 
one element per compatibility box. Any such grid of elements in the compatibility 
matrix we call a solution grid. The compatibility matrix encoding of G reduces 
DHC to search the encoding for solution grids. 

All elements of any solution grid equal 1, their indexes /i depend on their indexes 
i and their indexes v depend on their indexes j, onljB 



(1.4) 



1 | n = n(i), i/ = v(j), i,j = l,2,...,n} 



7 Here, we intentionally disregard the type of S and that particular value we fixed to it. 
^Index fi is a function of index i, and index u is a function of index j. 
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The following major properties of the elements' indexes follow from properties 11.31 
of the compatibility boxes: 



(1.5) 



* = J 
ii ^ 12 
h + h 



M(0 = "(j) 



- where all indexes are in their ranges. Due to these properties, any solution grid 
is a cyclical permutation of vertex indexes: 

( M (l) : /i(2) : M (3) : . . . : p(n)) = (i/(l) : i/(2) : i/(3) : . . . : 

In other words, it is a Hamiltonian cycle in G: 

M(l) -> M2) -> M3) -> ... ^ /i(n) -> 

Now, let's describe the guessing in the compatibility matrix's terms. It may be 
organized as follows: 

(1) Guess 7 is a box matrix with the same structure as C. In the matrix, all 
elements equal except those which create a possible solution grid in C, i.e. except 
those whose indexes create a set 

{(hOif*,") I a 1 = MO' v = Ki)) hi = M, ••-,«} 

- where functions and satisfy properties 11.51 All elements with these 
indexes equal 1. 

(2) There are n\ such guesses 7 in total. Let Y be set of all of them: 
(1-6) r = {71,72, ... ,7nl} 

(3) Checking conditions are the following relation: 

(1.7) 7 <C, 7 er 

Non-zero elements of any guess 7 which passed this test are a solution grid in C . 



1.3. Convex hull of all guesses. The following theorem clarifies the benefits of 
the compatibility matrix encoding. 

Theorem 1.1. Convex hull of auesses \1.6] is the following polytope: 

%ijfj,v — %jiv[i , •Eijiiv ^ 

- where i,j,fi,u= 1,2, ...,n, i ^ j, \i ^ v 

J2fj,= l, x ijfiv = Vjjvv 

- where i, j, v = 1, 2, . . . , n, i 7^ j 

En 
»=1, i^j ^ij'/if — Vjjvv 

- where j,fi, v = 1, 2, . . . , n, /i 7^ v 

^2u=i Vjjw — 1) > 

- where j = 1, 2, . . . , n 



(1.8) 
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Proof. System [l .81 may be described with the following box matrix of unknowns: 

/ 5*1,1 -X"l,2 • • ■ -X"l,ra ^ 
q _ -^2,1 ^2,2 ' ■ -^2,n 

\ Xn,l ■ ■ ■ Yn,n ) 

Diagonal boxes in B are the following nxn diagonal matrices: 

Yu = diag^j.i,!,! Vi,i,2,2 ■ ■ ■ Viiuu ■ ■ ■ Vi,i,n,n) 

Off-diagonal boxes in B are the following nxn matrices : 



Xi 



( 



V 



X 





•1- i 



J 



System 11.81 states the following relations between elements of 



i matrix B: 



(1) B is a symmetric matrix: Xji — Xfy Thus, all the below will be true in the 
horizontal direction as well as in the vertical direction; 

(2) For every fixed box column j and for every fixed column v in it: the total over 
/i = 1, 2, . . . , n, fi ^ of all elements Xij^ v does not depend on i. It equals to 



element 



of box F, 



33 i 



(3) For every fixed box column j, for every fixed column v in it, and for every fixed 
row fj, 7^ z/, the total over i ^ j of all elements Xij^ docs not depend on fi. It 
equals to element yjj VV of box Yjj; 

(4) For every box Yjj , its diagonal is a convex decomposition of 1 . 

Let's notice, system [TTH1 is consistent. For example, the following solution of the 
system we call a center: 



1 

n 



1 



n(n — 1) 

- where all indexes are in their ranges. The center minimizes Euclidean norm in 
polytope 11.81 Other obvious solutions of system 11.81 are guesses 11.61 The guesses 
maximize Euclidean norm in polytope 11.81 Thus, every guess is a vertex of the 
polytope. 

Now, let the above matrix B be a solution of system [l~8l It is easy to see that, 
for every fixed j cmd elements xij^i/^ i 7^ j and /i^c, constitute a (n — 1) x (n— 1) 
doubly stochastic matrix multiplied on yjjvv (the same is true for elements Xjiu^, 
i =/= j and fj, ^ v). Due to Birkhoff-von Neumann theorem [12] . that matrix is a 
convex combination of (n — 1) x (n — 1) permutation matrices multiplied on y 3JVV - 
Let's substitute these convex decompositions in B and write B as the appropriate 
sum. That sum is a convex decomposition of B over guesses [1.61 Thus, the guesses 
are the only vertices in polytope 11.81 □ 

1.4. The ATSP polytope. Due to theorem ll.il solution grids are those vertices 
of polytope 11.81 which satisfy inequality 11.71 But, that inequality can be expressed 
with the following equalities: 



(1.9) 
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- where indexes (io, jo, Mo, ^o) are indexes of all those elements in compatibility 
matrix C which equal 0. 

Theorem 1.2. Convex hull of all solution grids is volvtove \1.8l \1.9l 

Proof. Let's just notice that equations 11.91 cut from polytope 11.81 a chunk. That 
cutting is along the (io, jo, Mo, jy o)-th coordinate lines going through vertex points 
of the polytope, i.e. it does not create new vertices. □ 

Aggregated linear system 11.81 and 11.91 expresses the ATSP polytope. The sys- 
tem has polynomial size and can be solved in polynomial time. For example, 
Khachiyan's ellipsoid algorithm [Til H3] an d Karmarkar's interior-point algorithm 
[14j will solve this system in strongly polynomial time, because all its coefficients 
are or 1. 

Due to definition 11.21 of the compatibility boxes, system 11.91 explicitly involves 
an adjacency matrix of digraph G. Vertex relabeling of G will rotate chunk [L9l all 
over polytope 11.81 In other words, aggregated system fl .81 and Tl .91 is asymmetric in 
the sense of [TJ. Thus, theorem 11.21 may be seen as complimentary to Yannakakis' 
theorem [JJ: the ATSP polytope can be expressed by an asymmetric polynomial 
size linear program. 

Let's notice that, due to theorem ll.il any solution of aggregated system [L8l and 
II. 91 can be presented as a convex combination of guesses \l. 61 This decomposition is 
a P-problcm. Due to equalities 11.9} any guess in any such decomposition will be a 
solution grid, i.e. it will deliver a Hamiltonian cycle. Thus, theorem 11.21 efficiently 
solves DHC in both senses as a decision problem and as a search problem. 

Also, let's notice that the number of vertices in polytope 11.81 11.91 divided by n 
is the number of different Hamiltonian cycles in G. 

2. Asymmetric Traveling Salesman Problem 

Asymmetric Traveling Salesman Problem (ATSP) is a problem to find in any 
given weighted digraph a Hamiltonian cycle with minimal total weight. 

ATSP is a well known problem of combinatorial optimization [Il [T^[T6lfi7lfT^[T9l 
and many others]. Decision versiorQ of ATSP is a NP-complete problem [9]. Yet, 
theorem 11.21 allows an expression of ATSP by a polynomial size linear program. 

Let If be a weight function on given digraph G: 

W = (Wij)nxn ■ ^ Wij G (-oo, +oo], i,j = 1,2, ...,n 

- where Wij is the weight of arc from vertex i into vertex j (as usual, tuy = +oo if 
there is not any such arc). 

Theorem 2.1. The minimal total weight of Hamiltonian cycles in G is the solution 
of the following linear program: 

n 

(2.1) ^2 w lJ x l ^ v — > min 

- subject to constrains 1 1 . 8\ and 1 1 . 9[ 



9 The problem of existence of a Hamiltonian cycle with total weight in any given boundaries. 
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Proof. Due to theorem ll.21 G has Hamiltonian cycles iff program l2.1l has non-empty 
feasible set. Due to the same theorem, vertices of that set are Hamiltonian cycles 
in G, and the addends in criterion 12.11 pertain to the weights of those arcs which 
participate in the cycles. □ 

Linear program 12.11 has polynomial size and can be solved in polynomial time 
[ITl fT5l [2]. From the practical perspective, let's notice that we do not require 
weights Wij to be positive. 

Conclusion 

In this article, we reduced DHC to an asymmetric 0(n 4 )-size linear system and 
expressed ATSP by the appropriate linear program. 

The linearization was done by immersing 0(n 2 )-dimensional algebraic variety 
Oin R n . The immersion was done with the compatibility matrix encoding. 

Compatibility matrix is an encoding of DHC instances in the contradictions 
between relabeling options for vertex couples. The options can be computed in 
0(n 4 )-time with the brute force method. Analysis of the contradictions is a parallel 
testing of all guesses. The testing is the solution of linear system [L8l and [TT9l 

Let's notice that the role of those particular value and type of matrix S in 
inequality 11.11 was insignificant. So, our method can be directly applied to Subgraph 
Isomorphism problem (SubGI) in (multi) digraphs with (multi) loops |20j and to 
related optimizations. 

Because of the possibility of loops in SubGI (the diagonal elements in adjacency 
matrices Aq and S may be positive), the second property 1 1 . 31 has to be changed to 
inequality 

Ga < U n 

Strong inequalities in these relations cause additional constrains on variables ya^ 
in system 11.91 

That will be the only change for SubGI. 

There is a demo [21] . Using it, readers may try the reduction on their own 
examples. 
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